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ABSTRACT 

Project 1640 is a high-contrast imaging instrument recently commissioned at Palomar obser- 
vatory. A combination of a coronagraph with an integral field spectrograph (IFS), Project 1640 is 
designed to detect and characterize extrasolar planets, brown dwarfs, and circumstellar material 
orbiting nearby stars. In this paper, we present our data processing techniques for improving 
upon instrument raw sensitivity via the removal of quasi-static speckles. Our approach utilizes 
the chromatic image diversity provided by the IFS in combination with the locally-optimized 
combination of images (LOCI) algorithm to suppress the intensity of residual contaminating 
light in close angular proximity to target stars. We describe the Project 1640 speckle suppression 
pipeline (PSSP) and demonstrate the ability to detect companions with brightness comparable to 
and below that of initial speckle intensities using on-sky commissioning data. Our preliminary re- 
sults indicate that suppression factors of at least one order of magnitude are consistently possible, 
reaching 5cr contrast levels of 2.1 x 10^^ at 1" in the H-band in 20 minutes of on-source integration 
time when non-common-path errors are reasonably well-calibrated. These results suggest that 
near-infrared contrast levels of order « 10~^ at subarcsecond separations will soon be possible for 
Project 1640 and similarly designed instruments that receive a diffraction-limited beam corrected 
by adaptive optics (AO) systems employing deformable mirrors with high actuator-density. 

Subject headings: techniques: spectroscopic, high angular resolution - instrumentation: adaptive optics 
1. Introduction 
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Wavefront errors created by imperfections in 
the optical components (mirrors, lenses, beam 
splitters, ... etc.) of a high-contrast imaging in- 
strument manifest as a complex pattern of quasi- 
static intensity variations, or so-called "speckles" , 
in the focal plane at the science camera. Nomi- 
nally bright compared to the signal of a substel- 
lar companion or debris disk, speckle intensities 
must first be minimized by hardware, usually with 
an adaptive optics (AO) system, and then further 
suppressed through data processing. Since indi- 
vidual speckles have the same intrinsic width as 
the instrument point-spread- function (PSF), this 
latter step is crucial for distinguishing between 
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stellar artifacts and true companions. 

Speckle suppression is accomplished through 
differential imaging by exploiting differences be- 
tween the properties of residual starlight that 
reaches the detector and any incoherent radiation 
arriving from an off-axis source. To preserve com- 
panion flux, some form of "diversity" must be in- 
troduced into the system, whereby a sequence of 
images is recorded such that individual frames are 
trivially different from one another. Once mea- 
sured or calibrated, the image diversity is reversed 
using software to remove the speckles and reveal 
previously hidden companions. 

Image diversity can be achieved in a num- 
ber of ways. Examples include modulation 
of companion position with respect to speck- 
les through i nstrument rotation or field rotation 



(Krist 2007: Marois et al. 2008 



2010bl). m odu- 
lation of input polarizati o n state (Pott er 200S 



Qppenheimer et al.l l2008t iHinklev et al.l |2009), 



and mo dulation of th e actual target being ob- 
served ( Serabvn et al. 2010l : Crepp et al.l 2010l : 
Mawet et al.ll2009l ). The level to which residual 



starlight is removed depends on the modulation 
timescale compared to the speckle decorrelation 
timescale. Quasi-static speckle^ generally decor- 
relate in tens or hundreds of seconds as the result 
of changes in ambien t conditions and opti cal sys- 



tem alignment (e.g., iHinklev et al.l pOOTI )). Ide- 



ally, science images and PSF reference images are 
recorded simultaneously to minimize the influence 
of speckle pattern evolution. 

An alternative approach to discriminate be- 
tween companions and starlight is to generate 
diversity through the wavelength dependence of 
speckles. In this case, modulation of companion 
position with respect to speckles is achieved nat- 
urally through diffraction by acquiring images si- 
multaneously in multiple different spectral chan- 
nels. Chromatic differential imaging not only has 
the benefit of "freezing" the speckle pattern in 
time, but it also yields the spectrum of any com- 
panions, thus enabling study of their chemical 
composition and bulk physica l prop erties. First 
proposed bv ISparks and Ford (j2002l ). concurrent 
search and characterization may be accomplished 



^ Other families of speckles related to the atmosphere and 
adaptive optics system can deco rrelate on a much faster 
timescale l lMacintosh et al.ll2005f ). 



in practice using an integral-field spectrograph 
(IFS). 

Co mpared to other speckle suppression tech- 

nique s ( Absil and Maweti2009l : Qppenheimer and Hinklev 



20091) . chromatic differential imaging has a high 
duty-cycle efficiency and can help maximize in- 
strument scientific return by: accommodating a 
relatively wide bandpasH; allowing for the inter- 
changeable usage of science images and PSF refer- 
ences; maintaining an inner-working-angle that is 
independent of target coordinates relative to the 
observatory; avoiding smearing of the companion 
or dust disk PSF during an exposure; and obviat- 
ing the need to observe a nearby PSF calibration 
star with similar brightness and colors. Chromatic 
differential imaging may also be implemented in 
combination with other techniques. 

AO-assisted integral field spectroscopy has been 
used previo usly to study the comp anions orbiting 



mpa 

GO Lupi (iMcElwain et al.l l20o l AB Doradus 



(iThatte et al.ll2007l) . AB Pictoris (|Bonnefov et al 
20101 ). and recently the outer-planets of HR 



Bowler et all l2010fl . 



8799 (|Janson et al.l l2010at 
In this paper, we present the first on-sky ex- 
periments that combine an IFS with a coron- 



agraph (jHinklev et al. l2010bl) and speckle sup- 



pression using chrom atic differential imaging 
(jSparks and Fordll2002f ). To this, we also imple- 



ment the locally optimized combination of images 
(LOCI) algorithm which improves the signal-to - 
noise ratio of detections ()Lafreniere et al. 2007a ) . 

Project 1640 (hereafter, P1640) is a ground- 
based high-contrast imaging instrument that was 
recently installed and tested on the Hale 200- 
inch telescope at Palomar. The hardware incorpo- 
rates a near-infrared coronagraph and IFS and will 
soon receive a corrected be am from the PALM- 
3000 "extreme" AO svstem dHinklev et ahlboiObt 
Bouchez et al]l2009t ISoummeijl2005l ). P1640 has 
made several discoveries to date, including Alcor 
B and Zeta Virginis B ( Zimmerman et all 2010l : 
Hinklev et al.ll2010al) . Though faint compared to 
their host stars, these companions have masses of 
« O.25M0 and w 0.17Mq respectively, and were 
noticed in raw data. 



special case of chromatic differential imaging, called si- 
multaneous differential imaging (SDl), operates over a nar- 
row bandpass and has shown promise for high-contrast ap- 
plications involving me thanated sources as demonstrate d 
by the NICI campaign l|Biller et al.ll2010l : I5u et al.ll2010h . 
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In the following, we describe our technique 
for detecting companions having brightness and 
angular separation that places them beneath 
the noise floor set by speckles prior to data 
processing. This paper is accompanied by a 
companion paper by Pueyo et al. 2011 that 
discusses our method for subsequently extract- 
ing their spectrum. Combined, our results 
demonstrate the two principal utilities of using 
an IFS for high-contrast observations. These 
techniques are relevant to forth-coming instru- 
ments with similar designs, including the Gem- 
ini Planet Imager (iMacintosh et al.l l2006l ). VLT 
SPHERE ("Beuzit et al.' '2006^) , and Subaru Hi- 
CIAO ([McEl wain and SEEDS collaborationl2010[) . 
which will also use an IFS for chromatic differen- 
tial imaging and companion characterization. 

2. Data Pipeline Description 

2.1. IFS Kaw Data 

The P1640 IFS records 40,000 different spec- 
tra across the instrument field of view with every 
exposure. Each spectrum corresponds to an indi- 
vidual element of a 200 x 200 lenslet array and 
is mapped to a different location in the image 
plane. Each lens in the array has a pitch of sal9.2 
mas as projected onto the sky. The P1640 spec- 
tral resolution is 23 channels across the J and H 
bands, spanning wave lengths 1.10 < A < 1.75 /im 
(|Hinklevet al.ll2010bl) . \ law spectra are converted 
into a data cube consisting of two spatial coor- 
dinates and one color coordinate using the pro- 
cedure described in detail in Zimmerman et al. 
2010b (submitted). Each extracted cube is fiat- 
fielded, background-subtracted, and wavelength- 
calibrated based on observations of spectral stan- 
dards. The P164-0 speckle suppression pipeline 
(PSSP) starts with extracted and spectrally cal- 
ibrated data cubes. 

2.2. PSSP Code 

The PSSP is a custom program written in Mat- 
lab. This language was selected because it is flex- 
ible and contains a number of useful, fast, and ro- 
bust built-in functions. Moreover, Matlab allows 
for a trivial transition from running scripts in se- 
rial on one CPU core to running the same scripts 
in parallel on multiple cores, facilitating handling 
of the large number of images and calculations per 



image. To enhance reduction speed, we use the 
Bluedot super-computing cluster at the NASA Ex- 
oplanet Science Institute (NExScI) which incorpo- 
rates 16 nodes with 8 cores (2.5 GHz) per node. 
This feature permits pseudo real-time reduction 
while at the observatory, which enables follow-up 
observations of promising targets on subsequent 
nights within the same run to improve the fidelity 
of detections, spectra, and astrometry, and possi- 
bly take advantage of improved seeing conditions. 

2.3. Technique 

At its core, the PSSP involves precision stretch- 
ing and shrinking of individual images to utilize 
the spatial dependence of the speckle pattern on 
wavelength. This technique takes advantage of the 
fact that the position of a true companion is con- 
stant in each channel while speckles move radi- 
ally outwards from the star as a function of wave- 
length. For instance, the dependence of speckle 
position in the focal plane scales linearly with 
wavelength for phase aberrations created in the 
pupil plane. This form of diversity lends itself to 
differential imaging (SJT]). Stretching and shrinking 
images separated by several wavelength channels 
simultaneously aligns speckles and mis- aligns com- 
panions (Fig. 1). Upon subtraction, the highly 
correlated speckles are removed and companions 
survive, dramatically improving the effective de- 
tection sensitivity (SJ3]). 

The multitude of available reference images 
provided by the IFS naturally lends itself to appli- 
cation of the locally optimize d combination of iin - 
ages (LOCI) algorithm ( Lafr eniere et al.l l2007bl ) . 
Once images are optimally stretched for a given 
wavelength channel, we apply LOCI to construct 
an optimal reference, including both color and 
time, by weighting individual images with coeffi- 
cients based on a least-squares matching of local 
speckle intensities via matrix inversion. LOCI 
has been used with other speckle suppression 
techniques and provides higher signal-to-noise ra- 
tio detections compared to PSF su btraction in- 



volvi n g a single reference image ([Marois et al 
20081: iMawet et all l2009t iLafreniere et al.l l2009t 
Crepp et al.ll2010[) . Details of the algorithm are 



discussed below. 
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Fig. 1. — Occulted images of the star Alcor in individual IFS channels illustrating differential imaging with 
chromatic speckle diversity. Each image has a width of 3" x 3". (a) Single channel in the J-band after 
applying a high-pass Fourier filter but prior to image stretching. The position of Alcor B is marked with an 
arrow, (b) Same image radially stretched to optimally match the speckle pattern for a longer wavelength 
channel, (c) Single channel in the H-band. The speckle structure is highly correlated between the middle 
panels (see interior of white boxes), enabling precision PSF alignment and subtraction, (d) Differenced image 
displayed on a linear intensity scale. The signal of a real companion is incoherent with starlight and survives 
post-processing because its position is independent of wavelength. Non-differenced images are displayed on 
a square-root intensity scale for comparison of the speckle patterns. 



2.4. Program Outline 

An outline of the PSSP procedure is shown in 
Figure 2. We begin by retrieving extracted data 
cubes output by the Zimmerman et al. 2010b raw 
IFS conversion program. Data headers are then 
culled for pertinent information and the code is 
preconditioned by the user. This latter step in- 
volves specifying which images to use based on 
data quality; searching for astrometric fiducial 
spotfH; generating FOV, high-pass Fourier filter, 
and other alignment masks; setting basic parame- 
ters related to image cleaning and the LOCI algo- 
rithm; inserting artificial companions to assess the 
fidelity of the reduction; and initializing the num- 
ber of Matlab "workers" for parallel processing. 
Both occulted and unocculted images are loaded 
in order to perform relative photometry and cal- 
culate sensitivity in terms of contrast. 

Depending on target near-infrared brightness, 
several dozen cubes are generated for each star. 
This large data set, which we refer to as a hyper- 
cube when the time dimension is included (Fig. 

grid of wires placed in the pupil plane diffracts 
starlight into four bright spots centered on the 
star and may be used for astromet ric pur- 
poses llSivaramakrishnan and OppenheimeH |2006| : 

[Zimmerman et ani2oTol V 



3), consists of hundreds of images. In each channel 
of each cube, the spatial pixels are cleaned using 
standard techniques to correct for any remaining 
intensity outliers and negative values. Depending 
on whether we are searching for point sources or 
diffuse circumstellar material, like a dust disk, the 
images may be high-pass Fourier filtered to remove 
the low frequency component of the stellar halo 
and enhance the signal of both companions and 
speckles - the latter of which is useful for image 
alignment. The low frequency component of the 
image is removed using a quadratic attenuation 
profile, which we find has better noise properties 
than a linear or step-function filter. This stage of 
the procedure is time-consuming and so is sent to 
the NExScI Bluedot cluster for processing. 

Once cleaned, the images in each wavelength 
channel are stacked across the hyper-cube time di- 
mension to form composite images (Fig. 3). This 
step increases the companion signal-to-noise ra- 
tio in preparation for the LOCI algorithm. The 
speckles themselves are used to register images. 
We find that this approach is efficient and as reli- 
able as using the astrometric grid spots for align- 
ment, provided a custom gray-scale mask is used 
to prevent the coronagraphic spot or companion 
itself from biasing the result. A Fourier-based pro- 
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Fig. 2. — Block diagram of the P1640 speckle suppression pipeline (PSSP). Images are first extracted from 
raw IFS data via the procedure described in Zimmerman et al. 2010b to form cubes. The cubes are then 
processed using wavelength diversity and the LOCI algorithm to detect faint companions that are originally 
buried in speckle noise. Multiple CPU nodes from the NASA Exoplanet Science Institute (NExScI) Bluedot 
super-computer run custom Matlab scripts in parallel to enhance speed at several stages. 



gram registers the ima.ges wi th 0.01 pixel precision 
( Guizar-Sicairos et al. 2008fl . We find that allow- 
ing for small distortions in the position of speckles, 
by treating local image regions as a flexible mem- 
brane, results in only a marginal improvement in 
PSF matching. The number of temporal cubes 
used to form composite images depends on data 
quality and quantity. In general, cubes with suf- 
ficiently high speckle signal-to-noise ratio may be 
divided into a larger number of groups that are 
treated separately by the PSSP. The results may 
then be combined following image subtraction. 

The LOCI algorithm constructs an optimal 
PSF reference for each wavelength channel of the 
composite cube using neighboring frames in color 
and time. To generate sufficient wavelength diver- 
sity, and therefore minimize subtraction of com- 
panion light, the reference frames associated with 
a given composite image must be separated by 
several spectral channels, Ac = Co — Ci, where Co 
is the channel of the composite image with wave- 
length A = Ao and Ci is the channel common to 
a given set of reference frames. The minimum 
difference between channels is a function of com- 
posite image wavelength, Aq, spectral resolution, 
i5A, and field angle, 6*0, 

where a ~ 1 is the fraction of a diffraction width 
that PSF reference images are radially stretched 



I Ac, 



(1) 



and ^0 is expressed in the same units (unit- 
less). The optimal Ac range is determined by 
the competing effects of companion throughput 
and speckle correlation between channels. 

Fig. 4 displays the degree of correlation be- 
tween images as a function of wavelength and 
time following cleaning and optimal stretching and 
shrinking. All images are compared to a channel 
in the H-band corresponding to a wavelength of 
Aq = 1.55 /im in the first cube. A box is drawn 
around the reference images that may be sent to 
the LOCI algorithm for this case. The difference 
between two sets of wavelength channels is labeled. 
We generally include channels in the water-bands, 
located around 1-35 < A < 1.50 /zm, despite their 
comparatively lower signal due to atmospheric ab- 
sorption. It is likewise possible to use longer wave- 
lengths to achieve chromatic diversity. This figure 
demonstrates the utility of an IFS for speckle sup- 
pression. Speckles in stretched images are highly 
correlated in color - the degree to which is limited 
only by the intrinsic wavelength dependence of the 
instrument, but decorrelate after several minutes 
in time. 

Reference images that are closely spaced in time 
with Ac > Acjjjjjj are optimally stretched us- 
ing an iterative least-squares technique for speckle 
spatial matching that incorporates precision reg- 
istration via the Fourier-based method mentioned 
above. Images are resampled using bicubic inter- 
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Fig. 3. — Schematic of a P1640 hyper cube. Composite images shown to the left of the vertical dashed line 
are built by co-adding frames in each wavelength channel across the time domain to increase companion 
signal-to-noise ratio prior to performing PSF subtraction. In this example, images with red stripes are used 
to create the image in solid red. An optimal PSF reference image for each composite is constructed from 
the many available reference frames in the color and time domains (solid blue images). A small sub- frame 
is shown to indicate that PSF subtraction occurs locally. Composite images and usable reference frames are 
separated by Ac wavelength channels to minimize subtraction of companion signal. 



polation, where the intensity of an output pixel is 
calculated from a weighted average of the neigh- 
boring 4x4 pixels. A first guess for the optimal 
stretching parameter, S, is based on the wave- 
length of channels of interest using the scaling re- 
lationship, S = Xo/X, expected from speckle move- 
ments for phase aberrations located solely in the 
pupil plane. Here Xq again represents the wave- 
length of a non-reference (composite) image and 
A represents the wavelength of a reference image. 
From this starting point, we then adjust the scal- 
ing factor to converge to an optimal value accurate 
to AS « 10~^S. 

Figure 5 shows the optimal scaling factors, S, 
plotted as a function of wavelength for the same 
reference channel and dataset shown in Fig. 4. 
The theoretical relationship, S = Xq/X, is over- 
plotted for comparison. With each cube, we find 
that optimal scaling factors are systematically 
smaller than the expected Xq/X relation when A < 
Aq, but larger than Xq/X when A > Aq. This effect 
is understood intuitively by considering the radial 
displacement of speckles as a function of wave- 
length. The scaling relation, S = Xq/X, is derived 
for phase errors located in a pupil plane. In prac- 
tice, wave-front phase errors are also generated 
by out-of-pupil-plane optics and reflectivity and 
transmission non-uniformities which can trans- 



form into phase aberrations (jShaklan and Green 
20061 ) . Consequently, the radial position of a given 



speckle, 0{X), increments more slowly than the lin- 
ear dependence expected from first-order diffrac- 
tion theory and instead follows 9{X) < A 6*0 as a 
function of wavelength. 

Exceptions occur in the water-bands (A « 
1.4 /im) as a consequence of the cube extraction al- 
gorithm operating with low signal. In this regime, 
the spectrum in each spatial pixel becomes biased 
towards the red-end of the J-band and also blue- 
end of the H-band (see Zimmerman et al. 2010b 
for details). As a result, the cube extraction algo- 
rithm creates a spatial pattern that systematically 
mimics an adjacent spectral channel, thus creat- 
ing the bump seen in Fig. 5. This systematic 
error does not effect PSF subtraction, but it does 
artificially modify the spectrum of companions in 
the water-bands, requiring calibration. The rms 
value of S averaged across the spectral bands is 
0.003, comparable to image alignment precision. 
As a result of these effects, the PSSP recalculates 
the optimal scaling factor for each image in each 
cube for each star to maximize performance. 

Our implementation of the LOCI algorithm op- 
erates on each spatial pixel in turn. This step is 
likewise time-consuming and so is also sent to the 
Bluedot cluster. Local image coefficients are gen- 
erally assigned based on the flux inside of a box 
of width 5x5 PSF widths. Fig. 6 shows op- 
timization of the signal-to-noise ratio as a func- 
tion of the LOCI box width (referred to as the 
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Fig. 4. — Cross-correlation peaks in color and time for the star Alcor after optimal image stretching. The star 
is occulted by the coronagraph, the cubes have been high-pass Fourier filtered, and Alcor B has been masked. 
All images are compared to the Ao = 1.55 /im channel in the first cube. The legend shows the amount of 
time elapsed in minutes between exposures. Speckles are highly correlated as a function of wavelength (see 
also Fig. 1), thus justifying the use of processed IFS data for PSF reference subtraction. When constructing 
an optimal reference image via the LOCI algorithm, several channels must be skipped (given by Ac) to 
prevent subtraction of any companions. We nominally use cubes acquired within a certain interval of time. 
These images are identified as enclosed by a box. Different channels have different reference images. In this 
example, there are 18 highly correlated images available that provide sufficient wavelength diversity for the 
1.55 fim channel. More images may be used in general since the LOCI algorithm will assign their associated 
coefficients an appropriate (lower) weight. 



0-zone in Lafreniere et al. ( 2007bl )) for three ar- 
tificial companions inserted into the Alcor data 
set. Their angular separation from the star are 
0.40", 0.46", 0.52" respectively as shown in Fig. 7 
and discussed in more detail below. Once com- 
plete, the program co-adds the reduced data in the 
color dimension. The J and H-bands are separated 
at A = 1.40 /im to avoid unnecessary averaging in 
case a companion is particularly red or blue. 

3. On-Sky Demonstration 

P1640 has been used on several observing runs 
at Palomar in recent semesters. In this section, 
we present our preliminary on-sky results, demon- 
strating the ability to recover real and artificial 
companions that are faint compared to the inten- 
sity of quasi-static speckles. We characterize the 
signal-to-noise ratio of detections and show our 
sensitivity in terms of contrast. 



3.1. Alcor 

P1640 discovered an M-dwarf orbit ing the 

nearby, bright A-star Alcor in March 2009 (" Zimmerman et al 
i2010.) . With an angular separation of ?sl" and 
brightness ratio of 6 magnitudes, the w O.25M0 
companion was identified in stacked cubes prior 
to speckle suppression|f| In this section, we use 
the Alcor data set to demonstrate the ability to 
detect artificial companions that are fainter and lo- 
cated closer to the primary star where quasi-static 
speckles are bright as a result of non-common-path 
wave-front errors. 

To test the PSSP and also calculations of con- 
trast f iJ3.3p . we inserted artificial companions of 
various brightness and angular separation into the 
Alcor on-sky data. The companions were given 
PSF widths and relative intensities that match Al- 
cor B in each channel. To minimize PSF smearing 



^Extracted IFS data cubes may be converted into a movie 
that sequentially displays images in wavelength space, al- 
lowing for the visual identification of candidate companions 
whose brightness is comparable to the intensity of speckles. 
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Fig. 5. — Image scaling (stretching/shrinking) factors as a function of wavelength for all IFS channels 
calculated for Ao = 1.55 /^m using several different cubes of the star Alcor. If all aberrations where phase 
errors that occurred in the pupil plane, the optimal scaling factors would follow S = Aq/A, shown by the 
dashed line. Since aberrations are also generated in amplitude and by out-of-pupil-plane optics in phase, 
S < Ao/A when A < Aq, and 5 > Aq/A when A > Aq (see discussion). Image channels located between the J 
and H filters, near A « 1.35 — 1.50 /im, have a lower signal as they are affected by atmospheric water-bands, 
leading to a systematic offset and larger dispersion in the values by which images are scaled as a consequence 
of the cube extraction algorithm. We nominally use channels located in the water-bands to build PSF 
reference images for bright stars. Optimal scaling factors are therefore calculated on an individual basis. 
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Fig. 6. — Signal-to-noise ratio versus LOCI (outer) box width for three artificial companions injected into 
the Alcor data set shown in Fig. 7. The optimal LOCI box has a width of « 5 \/D, reaching a compro- 
mise between speckle suppression and companion throughput. The companions have equal brightness and 
separations of 0.40", 0.46", and 0.52". More wavelength diversity (Ac > 6) is required to detect companions 
with separation < 0.4". 
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Fig. 7. — Stacked H-band images of the star Alcor before (left) and after (right) running the PSSP. Three 
artificial companions have been injected into the on-sky data with separations of 0.40", 0.46", and 0.52". 
The companions are 1 magnitude fainter than Alcor B in each channel and detectable only after using 
the wavelength diversity provided by the IFS. The dipole pattern is evident for the inner-companions as a 
consequence of image stretching. Alcor B does not exhibit the dipole pattern because it is masked in PSF 
reference images to facilitate image registration. 



in color and time, each companion was injected 
into each image channel of each cube by using Al- 
cor B as an astrometric reference point. An opti- 
mal LOCI (outer) box size of 5 x 5 PSF widths was 
used to maximize the signal-to- noise ratio (Fig. 6). 
Fig. 7 shows H-band images of the star before and 
after application of the PSSP using three artifi- 
cial companions that are each 1 magnitude fainter 
than Alcor B with separations of 0.40", 0.46", and 
0.52" respectively. 

The three artificial companions are completely 
buried in speckle noise in pre-processed cubes and 
detectable only after using the wavelength diver- 
sity provided by the IFS. After application of 
PSSP, each display the tell-tale dipole pattern - 
a closely spaced positive and negative image - as 
a result of reference image stretching and subtrac- 
tionl3 Alcor B does not exhibit the dipole pattern 



^The inclusion of both image stretching and shrinking would 
yield a tri-polc pattern ( Sparks and Ford 2002). It is pos- 
sible to further build companion flux by subtracting com- 
posite images from reference images, and then rescaling and 
stacking the results. This too creates a tripole pattern, but 
it also smears the PSF and can degrade astrometric preci- 
sion. 



because it is masked in each PSF reference image 
to facilitate precision image registration. The arti- 
ficial companions have a final signal-to-noise ratio 
of 13, 24, and 31 for respective separations of 6* = 
0.40", 0.46", and 0.52" as shown Fig. 6. Alcor B 
is unambiguously recovered with a signal-to-noise 
ratio gain from 11.9 to 113.4. 

The artificial companions have a PSSP through- 
put of 42% (inner), 35% (middle), and 51% 
(outer), where throughput is defined as the ra- 
tio of the total companion fiux before and after 
data processing averaged over each wavelength 
channel in the H-band. Throughput is a func- 
tion of the local speckle noise and angular separa- 
tion and must be calibrated to accurately recover 
the spectro-photomctric signal of companions that 
are fainter than stellar artifacts prior to applying 
PSSP. We find that broadband aperture photom- 
etry measurements typically result in 10% errors 
following calibration, with larger uncertainties oc- 
curring closer to the star. More precise measure- 
ments may be obtained by selecting PSF reference 
images that minimize fiuctuations in the speckle 
subtraction residuals, but even 5% photometry is 



challenging given the many inherent biases present 
in LOCI data s ets tha t must be accounted for (e.g., 
Marois et al. I (l2ninal )'l. The topic of precise spec- 
tral extraction is discussed in detail in Pueyo et 
al. 2011. 

The bright ring surrounding the coronagraph 
occulting spot results from insufficient wavelength 
diversity as a consequence of the slower outward 
movement of speckles with spectral channel in the 
regions very close to the star. This signature may 
be used to qualitatively assess the quality of ref- 
erence image alignment. Residual starlight is dis- 
tributed symmetrically, but is more intense on one 
side of the coronagraph as a result of tip/tilt errors 
experienced in several data cubes. 

No real companions closer to the Alcor pri- 
mary were detected. We find that artificial com- 
panions injected with a yet fainter signal, by fac- 
tors of «2.6, 4.8, and 6.2 at separations of = 
0.40", 0.46", and 0.52", are recoverable respec- 
tively. This result is consistent with Fig. 6 and 
was also used to check PSSP contrast calculations 
(see §3.3|) . Alcor is a member of the Ursa Ma- 
joris moving gr oup which has a purp o rted age of 



400-6 00 Myrs (ICastellani et al.M2002l: King et al 



2003f) . Using the IZimmerman et al.l (|20ld ) H- 



band photometry of the primary and secondary 
along with our contrast measurements, we are able 
to rule out additional companions with ma ss m > 
77Mj and separation > 0.52" ((Baraffe et al 



20031) 



3.2. FU Orionis 

FU Orionis (Ori) is a prototypical young stel- 
lar object after which a class of enigmatic low- 
mass stars exhibiting signs of disk accretion is 
named. Such objects experience outbursts result- 
ing in 4-5 magnitudes of visual brightening on 
time scales of months with decay over decades 
sometimes punctuated by sm aller scale episodes 
( Hartmann and Kenvon 19961 ). FU Ori itself has 
a known stellar companion separated by 0.5" 
(|Wang et al.ll2004l ). It is thought that many FU 
Ori objects are members of binary or hie rarchical 
triple systems (jReipurth and AsDinll2004l) . 

We observed FU Ori on March 17, 2009 to bet- 
ter characterize the companion which we measure 
to be 5.2 magnitudes fainter than the primary in 
the J-band and 4.4 mags fainter in the H-band. 



Given the level of contamination, FU Ori B does 
not currently have a published J-band spectrum 
as it is mixed with speckle noise. Indeed, it was 
originally discovered us ing PSF subtract ion of a 
nearby calibration star (|Wang et al.ll2004l ). 

Fig. 8 shows combined JH images of FU Ori 
before and after running the PSSP algorithm. 
The companion is just perceptible in pre-processed 
data. Once applied, the PSSP unambiguously re- 
covers the companion as indicated by the bright 
positive and negative signature. Analysis of indi- 
vidual channels shows that we are able to improve 
from a 1.5(t to 16.4cr detection at A = 1.22 /im and 
2.3(7 to 22.9(7 detection at A = 1.58 /im, where a 
is the local standard deviation in speckle noise in- 
tensity. The J and H-band spectrum of FU Ori 
B will be presented in a subsequent paper. This 
experiment explicitly demonstrates the utility of 
using an IFS for speckle-suppression applications. 

3.3. Measured Contrast 

The PSSP provides sensitivity measurements 
for each star reduced. Contrast is calculated 
at several stages during the reduction, including 
with raw (pre-processed) data, high-pass Fourier- 
filtered data, and speckle suppressed data. Speckle 
noise is found by measuring the local standard de- 
viation of intensity variations occurring in a box of 
size 5x5 PSF widths, where flux levels are calcu- 
lated on the scale of a single PSF by averaging over 
several pixels. Speckle noise values are then com- 
pared to the stellar peak intensity measured in un- 
occulted frames to calculate the relative brightness 
of directly detectable companions. We also esti- 
mate the noise floor set by photon arrival statistics 
by measuring the square-root of the mean flux over 
the same regions. We nominally quote 5a contrast 
levels which correspond to a signal-to-noise ratio 
where companions become noticeable by eye0 

Resuhs for the V = 6.7 star HD 204277 
from data taken on June 28, 2009 are shown in 
Fig. 9. Local contrast levels using raw, Fourier- 
filtered, and speckle-suppressed images have been 
azimuthally averaged (median). We find that 



^Speckle noise follows a Rician spatial distribution 
and speckles with inten s ity up to RilOcr may occu r 
l lAime and Soummej [20041 : iFitzgerald and Grahaml |2006|) . 
An IFS however can help to discriminate between unusu- 
ally bright speckles and real companions by comparing the 
spectrum of the candidate to that of the star. 
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Fig. 8. — Stacked JH images of the star FU Ori before (left) and after (right) running the PSSP. The 
companion, FU ORI B, is located within the AO control box and has a brightness comparable to the 
scattered-light halo of the primary. Its unambiguous detection following post-processing demonstrates the 
principal of using an IFS to discriminate between speckles and faint off-axis sources. 



residual quasi-static wave-front phase errors of 
RilOO nm rms limit raw (pre-processed) sensitivity 
to 10^'^ at the coronagraph inner-working-angle. 
This result is consistent with similar measure- 
ments using the PHARO instrument prior to ap- 
plication of the modified Gerchberg-Saxt on phase 



retrieval algorithm ( Burruss et al. 



High- 
pass filtering the data removes the low-frequency 
component of the stellar halo (pedestal) resulting 
in a factor of several improvement, most notably 
at separations exterior to the AO control region 
at «0.6". Finally, chromatic differential imaging 
operating in tandem with the LOCI algorithm re- 
duces remaining noise by an order of magnitude 
closer to the star. In the case of HD 204277, this 
results in a 5 cr H-band contrast of 2.1 x 10"^ at 
1" with 20 minutes of on-source integration time. 

These are the deepest contrast levels yet 
achieved at Palomar in the H-band. Sensitivity is 
currently limited by calibration of non-common- 
path wave-front errors and instrument transmis- 
sion which ultimately governs the level of photon- 
noise and signal-to-noise ratio of speckles. We 
find that contrast continues to improve with ad- 
ditional exposures, but slower than a square-root 
relationship in time in regions very close to the 
star. 



4. Summary &: Concluding Remarks 

Central to the issue of generating contrast lev- 
els sufficient to directly image extrasolar planets 
is subtraction of (unavoidable) quasi-static speck- 
les that arise from instrument optical aberrations. 
Along with the ability to characterize companions, 
an integral field spectrograph naturally provides 
the capability for suppressing this dominant noise 
source. The recently commissioned Project 1640 
instrument at Palomar is the first to use an IFS 
for high-contrast imaging in combination with a 
coronagraph and adaptive optics. 

We have written a custom program, the Project 
1640 speckle suppression pipehne (PSSP), to 
improve instrument raw sensitivity via post- 
processing. The program utilizes a version of 
chromatic differential imaging, or spectral decon- 
volution, to perform point-spread-function sub- 
traction, as was first proposed by Sparks & Ford 
2002. To this, we also use the multitude of images 
provided by the IFS to construct PSF references 
via the locally opt imized combination of images 
(LOCI) algorithm (jLafreniere et al.ll20Q7bl) . 

We have quantified the degree of spatial corre- 
lation between speckle patterns in adjacent P1640 
hyper-cube images as a function of color and time. 
We find that spectral image diversity provides PSF 
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Fig. 9.— P1640 contrast levels in the H-band for the star HD 204277 taken on June 28, 2009 UT. The star 
is occulted by the apodized pupil Lyot coronagraph. PSSP consistently provides a gain of at least 1 order 
of magnitude in sensitivity at subarcsecond separations. 



references that match science images at a level 
comparable to, if not better than, observations of 
a nearby calibration star, depending on over-head 
time and system flexure. Analysis related to image 
scaling factors, throughput, and signal-to-noise ra- 
tio well-match theoretical expectations and results 
from other speckle diversity techniques that also 
use LOCI, once relevant parameters are optimized. 
Further, we are able to perform pseudo-real-time 
data processing while at the observatory using the 
Bluedot super-computing cluster at NExScI, 

We have quantified the improvement afforded 
by the PSSP with three different data sets. In- 
jecting faint companions, recovering known com- 
panions, and calculating contrast levels before and 
after performing speckle suppression each indicate 
that on-sky sensitivities consistently improve by 
at least one order of magnitude at angular sep- 
arations within the AO control region following 
post-processing. At larger separations, the in- 
tensity of residual starlight approaches the limit 
set by photon noise. Our sensitivity currently 
reaches 5cr contrast levels of « 2.1 x 10~^ at 1" in 
the H-band with 20 minutes of integration time. 
This result is comparable to other experiments at 
Palomar and is currently li mited by calibration 
of non-common-pa th errors ( Burruss et al. 2010l: 



Crepp et al.l |2010( ) . Preliminary tests using an 



interferometric wave-front calibration unit show 



promise to reduce these errors by an order of mag- 
nitude when operating in tandem with the PALM- 
3000 "extreme" AO system, suggesting that con- 
trast levels of order w 10^^ inside of an arcsec- 
ond are feasible with integr ation times of 1 - 2 hou rs 
(Shao et al. 2011, in prep., iBouchez et aD(|2009h ). 

IFS data provides the ability to play movies 
by sequentially stepping through each image in 
the color dimension. We have performed exper- 
iments carefully comparing our ability to detect 
faint companions using these color-movies and the 
PSSP. We find that movies perform remarkably 
well, and that the PSSP generally provides an ad- 
ditional factor of 2-3 in effective contrast in com- 
parison. This gain is most noticeable in close 
proximity to the star, just exterior to the corona- 
graphic spot, and is largely a result of implement- 
ing the LOCI algorithm to improve PSF match- 
ing. We play color-movies as part of our stan- 
dard reduction package to identify companions 
with brightness comparable to raw speckle noise. 
Stepping through wavelength channels is likewise 
helpful with differenced images. 

It is necessary to calibrate the effects of partial 
subtraction resulting from LOCI to accurately 
measur e the relative intensity of a faint com- 
panion jMarois et al. 2008; Thalmann et a l.ll2009t 



Currie et al...2010;. Bowler et al.,.2010;,Janson et al 
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2010bl: iMarois et al.ll2010a|) - This is nominally ac- 
complished by injecting artificial companions with 
known brightness into each data cube. Our pho- 
tometric analysis indicates that 10% uncertainties 
are typical for angular separations within the AO 
control region and that larger uncertainties, by 
factors of several, may occur very close to the 
star where only marginally sufficient wavelength 
diversity is achieved. 

High-contrast astrometry is likewise challeng- 
ing. The inclusion of astrometric grid spots aids 
with locating the position of the occulted star, but 
aggressive differential imaging techniques modify 
the shape and fiux distribution of the companion 
PSF, degrading precision. Our preliminary analy- 
sis using artificial companions indicates that sys- 
tematic errors of size w 0.5 spatial pixels are com- 
mon, even when using companion masks in ref- 
erence frames. Carefully chosen reference frames 
can help to mitigate these effects, but 1 mas as- 
trometry using near-Nyquist sampled data (~2-4 
spatial pixels per diffraction width in this case) 
may prove to be prohibitive unless more advanced 
algorithms are developed. 

This paper demonstrates one of the two impor- 
tant benefits provided by an IFS for high-contrast 
imaging: the automatic generation of useful PSF 
references for chromatic differential imaging and 
the detection of faint companions. A subsequent 
paper by Pueyo et al. 2011 will discuss the chal- 
lenges and possible solutions for accurately ex- 
tracting their spectra. 
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